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ABSTRACT 


Quantitative analytical procedures for relating selected water 
quality parameters to the characteristics of the backs cat tered 
signals, measured by remote sensors, require the solution of the 
radiative transport equation in turbid media. In this paper, we 
present an approximate closed form solution of this equation and 
based on this solution discuss the remote sensing of the sediments. 
The results are compared with other standard closed form solutions 
such as quasi-single scattering approximations. 
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1.0 INTRODUCTION 


The emergence of new environmental regulations restricting 
coastal zone developments combined with the desirability of developing 
coastal areas for energy and other important activities, and finally 
expansion of coastal jurisdiction to a two hundred mile zone, make 
fast response monitoring capabilities, possible by remote sensing 
applications, an attractive option. The question may be raised, why 
remote sensing is not employed on a larger than presently utilized 
scale. The answer lies simply in the inadequacy of quantitative 
procedures for obtaining water quality information from remotely 
sensed data. Presently several governmental, and industrial programs 
are addressing these inadequacies. Among these is a laboratory 
program at NASA-Langley Research Center (LaRC) . The basic experimental 
setup and the goals of this program are described below. 

The purpose of the LaRC laboratory program is Co determine 
whether the spectral characteristic of the upwelling radiance can be 
related to the amount and properties of the particulates and other 
pollutants. In the LaRC experimental setup, the beam of a solar 
simulator illuminates a large water tank filled with turbid water; 
the water turbidity in the tank is generated by the stepwise intro- 
duction of known amounts of sediments and other materials of known 
properties. An overhead detector system including a spectrometer, 
electronics, and a camera, measures the strengths and the characteristics 
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of the upwelling radiance signal in the wavelengths range of 300-1100 
nm. 

The Metrek Division of the Mitre Corporation has been assisting 
LaRC in the area of radiative transfer modeling. In the past, Metrek' s 
modeling approach has been based on Monte Carlo simulation techniques 
[|l,2]. These techniques' were originally selected to avoid mathe- 
matical complexities involved in the analytical solution to the radia- 
tive transport equation resulting from the inclusion of the geometry 
of the tank, in which the LaRC's laboratory experiment is conducted. 
Recent application of Metrek' s Monte Carlo simulation model has shown 
that the inclvision of the tank geometry may not be necessary for high 
turbidity levels (attenuation coefficient a 2: ^ meter of 'interest 

to the LaRC program. Other factors which may lead to the same 

/ 

conclusion with regard to the tank geometry considerations are: 

• Small reflectivity of the tank boundaries 
(~ 3.0, percent) 

• Large dimensions of the tank (diameter =2.6 meter, , 
height = 3 meter) 

• Small footprint of the overhead detector (circular 
spot of radius = 3.3 cm). 

The smallness of the footprint of the detector permits another 
simplifying situation. With this, it may be assumed that the incident 
radiation (spot size 30 .cm) behaves like a plane wave as far as the 
radiation emerging through the footprint of the detector is concerned. 
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The simplified geometry, and the plane- wave assumption, described 
above, allow the derivation of an approximate solution to the radia- 
tive transfer equation in turbid water. 

The present report discusses this closed form solution and 
illustrates the variations in the upwelling radiance for' situations 
which. are of interest in remote sensing applications. In order 
to establish common notation the optical parameters- of the turbid 
water are described- in the next section. 

1.1 Optical Parameters of Turbid Water 

In the absence of polarization the- following, parameters are 
necessary for optical characterization of turbid water: 

• Total absorption coefficient, a. 

• Total scattering coefficient, s. 

These coefficients have the dimension of meter The ;attenuation 
coefficient, a, is the sum of the absorption and scattering .coeffi- 

cients. The single scattering .albedo, is defined as the ratio of 

total scattering to the total attenuation coefficient. 

Another parameter' of' interest is the scattering phase function, 
P(0). This function specifies the angular pattern of the scattering 
of a collimated beam from an infinitesimal volume of turbid water. 

The scattering probability fimction P^(vi), for polar angle 

(6 = cos ^ y) is defined in terms of the scattering phase function 

by: 
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P(y) dy 


P(y) dy 


for all y between rl and +1. 

1..2 Summary 

The purpose of this work is to provide a closed form analytical 
expression which will enable the Langley Researchers to verify their 
experimental results, and/or to optimize their experimental para- 
meters without considerable theoretical calculations. 

In Section 2, the; mathematical, formalism applicable to the 
radiative transport process is established for the geometry of the 
light scattering relevant to the LaRC laboratory set up. experiment. 
The analytical solution to the equation of radiative transfer is 
then .derived for the backscattered radiance, as a function, of m , 

incident radiance and an appropriate representation of. the . 

■k 

scattering phase function P(y) in the form P(y) = F 6(y-l) + B. 

This representation- adequately displays the properties of the 
scattering phase functions for turbid, water (i.e.,.a strong peak in . 
the forward direction), and at the same time allows derivation of . 
a closed form solution to the radiative transfer equation. 

t . • 

. In Section 3 the results, on measured as well as calculated 
scattering phase functions have been summarized. The measured phase 


F, and B are constants , and 6 is the Dirac Delta' function. 


1-4 



functions represent a wide range of turbid water conditions. The 
calculated scattering phase functions are arrived at utilizing the 
Mie formalism for polydispersed suspensions- and the size distribution 
measurements provided through the LaRC Laboratory program. Acceptable 
values for B, are then obtained using measured as well as calculated 
phase functions results. It is found that B varies between 0.007 
(for a Feldspar particulate distribution calculated from the Mie 
formalism) and 0.0585 (for measured San Diego Harbor Water). F can 
be shown to be related to B by F = 2-2B (see Section 2.2). 

Section 4 presents the numerical results determined from the 
solution of the radiative transfer equation obtained in Section 2.0. 

The backscattered radiance is calculated for various values of lo 

0 

and B. “q the single scattering albedo and varies with wavelength-, 

A. Typical values of the a/s ratio at 550 nm may vary between 0.06 
(for quartz particles) and 0.20 for natural waters [ 3 ]. The a/s 
ratio, for sediments as a fiinction of wavelength may vary between 
0.20 at A = 550 nm and 0.50 at A = 370 nm [4) . The results generated 
in this work are compared with those of Gordon [s]., who used a quasi- 
single scattering approximation solution. For smaller values of 
B(s= 0.02), and m^(=; 0.80) the difference in the two results is of 
the order of 15 percent. At higher values of 0.95) and B(~ 0.06), 

the difference in the two results is significant. 

^Specifically lio = .... = (1 + -^4^) where A is the wavelength, 

o s(A) -f a(A) s(A)'^ 
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2.0 SOLUTION TO THE RADIATIVE TRANSFER EQUATION 

2‘.'1 Radiative Transfer Equation For An Idealized Geometry .Case 
The radiative transfer equation in any medium may be written as 

(6.) 


n . VI(r, n) +’ a(r) ,I(t,ri) = 


/ 


-dS2(n ) P(n,n ) i(r,ii ) 


where ‘I (r,n)-> (see Figure 2,1) represents the radiance of a given 

* 

-frequency, ,at the position .defined -by r and in the .direction desi:g- 
nated- by the unit vector ti. P(n,rt ) is the .scattering phase function. 


and n -is the scattering .angle,. P(n,ri ) does not .include any polari- 
zation 'effects and satisfies the condition: 


( 1 ') 


_1 

4n 



P(n,Ti ) 


A 

dfi(q ) 


1 


(la) 


.dJ2(n ) -denotes the differential solid angle about q . A' rigorous 
derivation of Equation (1) together with a discussion of its range" 
of validity may be found in the literature [7,8] . 

Equation (1) may be specialized for' an ideal case of ah incident 
plane wave. .in which the incident light is unpolarized and is normally 
falling on a homogeneous half space. 'The homogeneous half space is 
represented by an attenuation coefficient, a, -and a scattering ■ 
coefficient, s. Under" these conditions Equation (1) becomes.: 

- f / I(..r) .r (2) 
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U is the cosine of angle 0 shown in Figure 2-1. Extensions to 
the geometries other than the one considered here are straight 
forward and will not be discussed in this report, 

s 

Using the transformation az=T, and defining , we obtain 

m r 

(P + 1) / PCPjP") dy" . (3) 

2, 2 The Scattering Phase Function P(PiP'^) 

The phase function P(y,y'') for a variety of turbid waters have 
been studied in a previous report [l] . In general, these functions 
show a very strong peak in the forward direction. As a result of 
this the standard method of solution to the equation of radiative 
transfer, which is based on the expansion of P(y,y'') in terms of 
Legendre polynomials is not practical. Legendre polynomial expan- 
stions involve large number of terms and result in a large number of 
coupled differential equations. 

A representation which displays the properties of the phase 
functions for turbid water and at the same time enables us to 
seek an analytical solution to the radiative transport equation is 
given by the distribution: 

P(y,y") = F dCy-y-) + B (4) 

where S is the Dirac delta function, B < < F, and 

F + 2B = 2 (4a) 
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Equation (4a) is the result of substituting Equation (4) in 
Equation (la) . 

2.3 Treatment of Multiple Scattering 

Upon substitution of Equation (4) in Equation (3) it may be shown 

that 


91 . T/ 
y 9:^+ I(x 


0) F 

,y) 1 1- 


u) B 
0 


I(t,mO d]i' 


On taking the Laplace transform of both sides, we obtain 


(5a) 




t -| / ( 1 ) F \ (0 i5 

s - 1(0, y)J + f 1 - ^ )<!>(s,y) = <f.^(s) 


(0 B 


(5b) 


00 

here $( 3 ,^) = J I(x,y) e^® dx , and 


^>^(s)_ = J (j)(s,y) dy 
-1 


Rearranging terms in Equation (5b) results in 


(j)(s,y) = 


‘^0- 4>_(s) + yl(o,y) 
2 ^ 

1 + vs - {0 .F/2 

0 


(6a)' 


Integrating both sides w.r.,t. y we get 
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1 


1 M ' I ^0 , ] 


,dy 

» 

'b 

“ ^ 1 / / 

-r, v” 

Bu 

f “S’ \ f ( 

0) T \ 

0 

( 0 1 / I - 


1-2S 

[h - 2 + ":) / h ■ 

- 2 - 


Equation (6b) can be further simplified if it is considered that 
I(o,y) = 6'(y-l) for y > o 


(6b) 


= r(o,.y) to be. determined for y. < o 

This means that the radiance just below the water surface and traveling 
in the water is the same as the incident radiance. I(o,y) for y < o 
is determined from the relationship between and lXo,y) obtained 


from Equation (6a) by substituting s = - 
written as 


uj B r / ml 


Equation (6b) then reduces to : 





and may be 


( 7 ) 
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2.4 Quasi-Single Scattering Approximation 


Before we discuss a more complete solution ,to Equation (8) we 
show how quasi-single scattering approximation can be derived from 
the present formulation, and study the assumptions underlying this 
app roximati on . 

The quasi-single scattering approximation as developed by 
Gordon [5] corresponds to neglecting terms which include B in 
Equation 08) , l.e. , 


^ (s) = I 

o 


/ 


m' F 

II - 0 + s/ 


(9) 


which upon making use of Equation (7) results in: 

I 


T,''' 

Bu) 


I(o,y) = 


2(1 ‘ (1+1^1) 


( 10 ) 


A further assumption in' the derivation of the quasi-single scattering 
approximation is the substitution of B in Equation (10) by P(- | P j) • 
Considering . this the resultant quasi-single scattering approximation 
becomes 


1(0, y) 


m^P(-|y|) 

2(1 - 



( 11 ) 


The Fresnel transmission coefficient which takes into account the 
reflraction due to the transmission from water into the air has been 
left out of Equation (11) . 
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2.5 Derivation of a Simplified Multiple Scattering Solution 
Here the Equation in point is (8). The denotninator of Equation 
(8) displays complex singularity characteristics (because of the 
branches of the log- function) in the physical range of the variables 
of interest. The derivation of the solution of interest is thus 
based on .an analytical continuation technique. In what follows we 
discuss the characteristics of the various terms in Equation (8) , 


which may be re-written as: 



The observations which can be made at this point are: 


(1) G(s) is regular for Re(s) < 1 - 


0) F 
o 


(2) is regular for Re(s) > o (by definition). 


\ 

(3) [1 ^ +s ) has a zero at s 


1 - 


m F) 
0 


( 12 ) 
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The question to be answered now, concerns the zeros of the 

denominator of Equation (12). Two cases are to be distinguished (a) 
1 


u 


- 11 < B which is satisfied only by a scattering medium with 


negligible absorption, and (b) Sf— 1| > B.. Only case (b) has 


physical significance for the case of radiative transfer in turbid 
water and is discussed in this report. 


B(i3 

Define x(s) = 1 - -r— ^ In .1 

2s 


which may Be re-written as 




(15a) 


t(s) = 


1 - 


T, 

0) F 
0 


1 - 


2s 


In 


as 



0) F 

1 - 

^ H- 1“ 

2 ^ ® 


'O 

0) F 

^ T _ 

0 


+ 


1 - 0 ) 
c 

2s 


In 


1 - 

'C. 

0) F 

° + c; 

2 ^ ® 




0) F 

M _ 

0 

— c 


(15b) 


For s-»-o, the term in the capital brackets goes to zero like s , 


and the second term approaches the value 


1 - 0 ) 


1 - CO F/2 
o 


. This in 5 >lies that 
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T_(s) 


G(s) 


(18) 


The function. is regular for Re(s) > o, and t_(s) is regular and 

non zero for Re(s) > -3* Thus the left hand side of Equation (18) 
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is regular for Re(s) > o. Furthermore, G(s) is regular for Re(s) < 


I 

II T_|_(s) is regular and non zero for Re(s) < 3. Hence 

the right side of Equation (18) is regular for Re(s) < g. Since 
3 > 0 these two domains of regularity overlap, and ^in the strip 
o < Re(s) < g the right and left sides of Equation (18) are equal 
(see Figure 2.2). Therefore they are equal everywhere. From the 
principles of . analytical continuation each side of Equation (18) is 
equal to a function H(s) which is regular in the entire complex 
s-plane. Since G(s)/x_^(s) is bounded for |s[ then from 

Liouville’s theorem [lo] , the function H(s) = constant and 
Equation (18) results in. 




C 


0 


T_(s) 

1 “ (0 F/2 + s 
0 


(19) 


By letting s o in Equations (19) and (6b) we have 



( 20 ) 


( 21 ) 


Substitution of Equations (20) and (21) in Equation (19) results in 
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FIGURE 2.2 

DOMAINS OF REGULARITY 



or 



-p-i“ 


To put In T^(s), Re(s) > o in a form more convenient than in Equation 
(23) i we make use of the Canchy Goursat theorem (11). From the 
application of this theorem the path is shifted from the line -3 + iy 
(_ 00 < y < oo) to the imaginary axis indented as shown in Figure 2.3. 
The contribution from the semicircle at z=o goes to zero when sj^o. 
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1" 



,-g-ioo -ico 


FIGURE 2.3 

AN EXAMPLE OF A CONTOUR INTEGRAL 
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Since s^o in Equation (23) , therefore we obtain 



Let z = ix 



—CO 
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and 


T_(s) = Exp 


'Xj 

(1) F 

s(l 1-) 

TT 



In 


Bto 


1 - 


W 


1 " w F/2 tan W 
0 


w F'^ 
o 


dW 


C27) 


2 2 2 
Sin W + s Cos W 


To evaluate lim t (s) , -.we' first consider lim 

s o ■ s o 

from Equation (-24) is given by 




which 


In 


lim In T (s) = 


Bu 


1 - u F/2 + z 
1 - ^ In ' - 


2z 


1 - u)^F/2 - z/J 


S -5“ O 


2iri 


dz (28) 


-x« 


To take into account the contribution from the semicircle at z=o. 
Equation (28) may. be rewritten .as 



-ie 
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Let' z=-ix in the first integral, and z=ix in the third integral of 


the right hand side of Equation (29). It will result in 



Upon substitution of 2 = e 



we, obtain 



( 32 ) 


2^16 



Thus 


lim 
s o 




Bco 


1 - 


1 - 


'\j 

0) F 
0 

2 J 


-U 


1 - (0 F/2 
0 


\ — 


( 33 ) 


We may now return to Equation '(22), and substitute in this equation 
Equation (33). This results in 


<{»o(s) 


|1 - (0 F/2 
o 


I _ , 11 

\ °/J 1 - u F/ 


(34) 


F/2 + s 


Recall that (from Equation 7), I(o,ji) = 


'\j 

to B 
o 




1 - 


2p o ^ y 

all y between -1 and 0. This making use of Equation (34) gives us 


, for 


'X; / 

(i) B / 

1(0, y) = I 1 - S F/ 2 I ‘ (l - m 

2 O 1 0/1 0 


-3/2 T - 


1 - (0 F/2 
o 


(35a) 


1 - u 


1 - m^F/2' 

The function' r |- | may be evaluated from Equation (27) and 


results in 


T - 


1 - to F/2) 
o ’ 


= Exp 


Bto 


InJ 1 - 


W 


1 ^ W 

1 - (0 F/2 
o 


. 

. 


2 2 2 2 
VI Sin W + s Cos W 


'] 


dW (35b) 
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Finally, I(o,y) is given by 


u B / u F , 

I(o,u) = 1 --f-1 


-3/2 


1 - y 


In 


Exp 


ii 

IT 


1 - 


B(jU W Cot w 
0 

1 - w F/2 
o 


2 2 2 2 
y Sin W + s Cos W J 


(36) 


for all y between -1 and 0, The I(o,y) represents the backs cattered 
radiance received by the overhead detector in a remote sensing appli- 
cation. The integral involved in Equation (36) may easily be evaluated 
for a specific scattering angle. 

Upon examination of Equation (36), it is evident that for a given 
scattering angle y, the backscattered radiance depends on the scattering 

'\i 

phase function (P(y) = f(B,F)), and the single scattering albedo o)^. 

r\j 

The functi onal relationship between I(o,y) and in our result 
(Equation 36) is different than the functional dependence of I(o,y) 
on (jj^ in the restilts obtained from quasi-single scattering approxL- 
mation |_5J . The difference in the two results becomes large as 
takes the values close to 1. The numerical results obtained from 
Equation (36) and their differences from quasi-single scattering 
approximation are reported in Section 4. 
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3.0 SCATTERING PHASE FUNCTIONS 

In this section we have summarized (1) the available information 
on the measurements of the scattering phase functions, and (2) the 
calculated scattering phase functions utilizing the Mie formalism 
for polydispersed suspensions on the basis of size distribution 
measurements provided to us through the LaRC laboratory program. 

The compilation of this information is primarily done to determine 
the constants B, and F for the phase function F(ij). These constants 
are required to evaluate Equation (36) in Section 2. 

3.1 Measured Scattering Probability Functions 

Several in-site, as well as in vitro measured scattering phase 
functions covering turbid to clear water conditions have earlier 
been reported in Reference [2j . Reference ^2^ presents (1) observa- 
tions made in lake water [12], coastal waters Atlantic surface 

water [l4] , Pacific near coastal water [isj , Mediterranean [lb ] , and 
Sargasso Sea water [l?] , (2) measurements taken by the Scripps 
Institution of Oceanography [l8]j in deep clear oceanic water (Tongue 
of the Ocean)-, near shore ocean water (off shore of Southern 
California), and very turbid harbor water (San Diego .Harbor,) (3) the 
observation taken in vitro by Petzold [iS^ to determine the effect of 
adding scattering and absorbing materials in the water. 

From these observations, a set of scattering probability function 

P (0), have been obtained, and are shown collectively in Figure 3-1. 
h 
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Scattering Angle - 0® 


FIGURE 3.1 

SCATTERING PROBABILITY FUNCTIONS FOR NATURAL OCEAN WATERS 
AND FOR FRESH WATER THAT HAS BEEN FILTERED 
AND ARTIFICIALLY MODIFIED 




3.2 Calculated Scattering Probability Functions 

The extension of Mie theory to the case of polydispersed 
suspensions has been thoroughly discussed in Reference [2jj along 
with the computational methods used to calculate the scattering 
phase functions. Empirical cumulative size distribution for the 2 
samples, Ball Clay, and Feldspar were provided by NASA/LaRC. These 
two distributions differed radically from each other, and thus 
allowed a reasonable investigation of the effect of different size 
distributions on the scattering phase functions. For both the 
distributions, the index of refraction (M) was chosen to have real 
as- well as imaginary part. In order to determine the effect of 
absorption, the imaginary part of the index of refraction was set 
to zero. 

The scattering probability fimctlons obtained from Mie theory 
calculations utilizing these two distributions are presented in 
Figures 3-2 to 3-5. Two cut off limits, 100 pm and 10 pm were 
considered to be appropriate for the computations. The different 
cut off limits allow to investigate the effect of sediment settling 
in the water. 

3. 3 The Constants , B and F 

From the definition of the scattering probability function, 

Pj.(e), the scattering phase function F(p) is given by [see Section l.l] 
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. FIGURE 3.2 

VOLUME SCATTERING DISTRIBUTION FUNCTIONS FOR FELDSPAR 
(lOOjUlVi CUTOFF X = 500NIVI) 
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FIGURE 3.3 

VOLUME SCATTERING DISTRIBUTION FUNCTIONS FOR FELDSPAR 
(lOjLlM CUTOFF X=500NM) 
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J P(ji) dp 


^ 


, ■ p = cos ^0, and ^l<p<l. 


/ P(y) dp 
-1 


Substitution of P(p) = F6(p-1) + B (see Equation 4) results in 

1 


Pp(G) = 


J [F6(p-1) + b] du 
p _ F+B - Bp 

F + 2B F + 2B 


(37) 


The 


1 

condition, % J |^F6(p-l) + B j dp = 1 


implies 


-1 

F + 2B =2 


and 


B = 


2-2Fp(0) 

1+ p 


(38) 


The numerical values of B as determined from Pt,( 6) illustrated 

’ r 

in Figures 3-1 to 3-5 are presented in Tables 3-1 and 3-II.' Four 
water qualities, San Diego Harbor, Tongue of the Ocean, offshore 
Southern California, and fresh water + scattering + absorbing 
material x<rere chosen to calculate the range of the values for B in 
Table 3-1- San Diego Harbor water represents the most turbid water 
of all, and gives an upper bound on the value of B(»0.06) for the 
measured scattering probability functions. The upper bound on B as 
determined from the calculated scattering probability functions 
utilizing Mie theory formulation is 0.041. 
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TABLE 3-1 


B THE fraction BACKSCATTEKED DETERMINED FROM IN SITU AND IN VITRO EXPERIMENTS 



WATER QUALITY 

Angle of 
Incidence 
0 

y = Cos0 

' 

San 

Diego 

Harbor 

Tongue of 
The 
Ocean 

Offshore 

Southern 

California 

Fresh Water 
+ scattering 
+ absorbing 
material 

45° 

0.707 

0.0585 

0.0585 • 

0.0585 

0.0468 

60° 

0.5 

0.05333 

0.05333 

0.05333 

0.0266 

o 

o 

0.324 

0.0447 

0.0447 

0.0447 

0,0223 

o 

o 

CO 

0.173648 

0.0341 

0.0341 

0 . 0341 

0.0170 

90° 

0.0 

0.03 

0.03 

0.03 

0.010 
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TABLE 3-II 


B THE FRACTION BACKSCATTERED CALCULATED FROM MIE THEORY 


Sediments 



Feldspar 

Ball Clay 

Angle of 
Incidence 
0 


j 

Cos 0 

100 y 
cut off 

10 y 
cut off 

100 y 
cut off 

10 y 
cut off 

45° 


0.707 

.041 


.039 

.041 

..039 

o 

o 


O'. 50' 

.038 


.037 

.038 

.037 

O 

o 


0.342 

.03. 


.028' 

.031 

.029 

O 

o 

00 


0.173648 

CM 

O 


.017 

.021 

.018 

o 

o 


0.0 

.008 


.007 

.008 1 

.007 


4.0 RESULTS 


In this section we describe our numerical results on the back- 
scattered radiance determined from the solution of the radiative 


transport equation. Before this is done, however, we present the 
closed form analytical solution derived earlier and describe 
some of the parameters involved. The normalized backscattered 
radiance may be written as (see Equation 36): 



s s 

0 ) = Single scattering albedo = — = — ; — 

0 ® a s+a 


S = Total scattering coefficient 

a = Total attenuation coefficient 


a = Total absorbing coefficient 
p = Cos0, 0 is the scattering angle 
B, and F are the parameters used to represent the 
phase function P(p) = F5(p-1) + B 
The range of the values of B have been determined from the 
measured as well as calculated phase functions and is fotind to 
vary between 0.007 and 0.06. F is computed' from B by a 
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relationship F = 2-2B, determined in Section 2, 

The integration involved in Equation (39) was performed utilizing 
a computer code given in Appendix A. The computer program makes use of 
a Lobatto [ 19 ] integration algorithim. 

To investigate the effect of wavelength, the single scattering 

<\j 

albedo, u xras varied between 0.6 and 0.95. These limits on u were 
o o 

arrived at from the limits on the a/s ratio which is a function of 
the wavelength. The a/s ratio for sediments varies between 0.20 at 
A = 550 nm and 0.50 at A = 370 nm [4]. . Typical values of a/s ratio 
at 550 nm are between 0.06 (for quartz particles) and 0.20 for 
environmental waters . 

The results of the calculations for backscattered radiance per 
unit incident radiance are presented in Figures 4-lA to 4-lG, and 
indicate that the backscattered radiance received by an overhead 
detector is very sensitive to the changes in (1) the scattering phase 
function P(ij) (P(]i) = f(B,F)), and (2) the single ‘scattering albedo, 

•V 

(0 . 

0 

Another observation which can be made from Figure 4-lA to 4-lG 

9I(o,p) 3I(o,vi) . ^ . ..... 1 

xs that — — >> — ^ — - — , x.e. , for a gxven scatterxng phase 

o 

fimction, the incremental change in the backscattered radiance w.r.t. 

is much larger than the incremental change in the backscattered 

radiance w.r.t. the scattering angle, y. At the same time the function 
( 

31(0 ,y ^ varies slowly for smaller values of cu (up to m = 0.80), and 
3y 0 0 
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attains a large value for close to unity. For example, for 

= 0.95, and P(p) = 1.88 6(p-l) + 0.060, the received backscattered 
radiance is as large as 80% of the incident radiance. 

4.1 Comparison of Metrek's Results with Quasi-Single Scatterin 


Approximation 


The quasi-single scattering approximation results were earlier 
reported by Gordon [5]. The backscattered radiance from his results 
is given by (see Equation 11) 


lQ(o,ii) = 


CO P(-- u ) 




As described in Section 2.4, the above equation is equivalent to 


Iq(o,p) = 


2 1 - 


1 - p 


-l<p<o 


In that case, from equations (39) and (40), we obtain 
I(o,.p) = FACTOR • Ip(o,p) 


where 


m^F 3/2 


FACTOR ^ 


IT, 

t/ 


Bco W Cot W 
In - 1 ° ^ . 

1 j 1 

2 ^2 2 ~ ^^ 2 ) 

U Sin W + S Cos W J 


The FACTOR was evaluated for'co^ ranging from 0^60 to 0.95, and B 
ranging from 0.007 to 0,06. The results are presented in Tables 4-1 to 4-VIl. 
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TABLE 4-1 


FACTOR FOR B = 0.007 AND F/2 = 0.993 


\o 
U X 

0.60 

0.65 

0.70 

0.75 

0.80 

0.85 

0.90 

0.95 

0.1 

1-017 

1.021 

1 .027 

1.034 

1.046 

1.065 

1 . 104 

1:224 

0.2 

1.018 

1.022 

1.028 

1-035 

1.047 

1.067 

1 - 108 

1.233 

0-3 

1.018 

1.022 

1.028 

1.036 

1.049 

1.069 

1.111 

1.239 

0.4 

1.018 

1.023 

1.029 

1.037 

1.049 

1.070 

1.113 

1.244 

0.5 

1.019 

1.023 

1.029 

1.038 

1.050 

1.071 

1 . 114 

1,247 

0.6 

1.019 

i .023 

1-029 

1.038 

1.051 

1.072 

1.116 

1.250 

0.7 

1.019 

1.024 

1.030 

1.038 

1.051 

1.073 

1.117 

1.253 

0.8 

1.019 

1.024 

1 .030 

1-039 

1.052 

1.073 

1 . 118 

1.255 

0.9 

1.019 

1.024 

1.030 

1.039 

1.052 

1.074 

1.119 

1.257 

1.0 

1.019 

1.024 

1.030 

1.039 

1.052 

1.074 

1.119 

1.259 


31 - 


TABLE 4-II 


FACTOR FOR B = 0.010 AND F/2 = 0.990 


Vo 

V \ 

0.60 

0.65 

0.70 

0.75 

0.80 

0.85 

0.90 

0.95 

0.1 

1.024 

1.030 

1.038 

1.049 

1.066 

1.093 

1.150 

1.324 

0.2 

1.025 

1.031 

1.039 

1.051 

1.068 

1.097 

1 . 155 

1.338 

0.3 

1.026 

1.032 

1.040 

1.052 

1.070 

1.099 

1 . 159 

1.347 

0.4 

1.026 

1.033 

1.041 

1-053 

1,071 

1.101 

1 . 162 

1.354 

0.5 

1.027 

1.033 

1.042 

1.054 

1.072 

1.103 

1.165 

1.360 

0.6 

1.027 

1.034 

1.042 

1.054 

1.073 

1,104 

1.167 

1.365 

0.7 

1.027 

1.034 

1.043 

1.055 

1.074 

1.105 

1.169 

i .369 

0.8 

1.028 

1.034 

1.043 

1.055 

1.074 

1.106 

1 .170 

1.372 

0.9 

1.028 

1.034 

1,043 

1.056 

1.075 

1.106 

1.171 

1.375 

1.0 

1.028 

1.035 

1.043 

1.056 

1.075 

1.107 

1.172 

1.377 
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TABLE 4 -III 


'© O' 

.taj ^ 

S3 

<Q 

Sff 

H O 
S tS? 


a> 


FACTOR FOR B = 0.020 AND F /2 = 0.980 


\ 0 ) 

N. 0 

li \ 

0.60 

0.65 

0.70 

0.75 

0.80 . 

0.85 

. 0.90 

0.95 

Oil 

1.049 

1.061 

1.077 

1.099 

1.133 

1-190 

1.306 

1.680 

0.2 

1.051 

1.063 

1-080 

1.103 

1.138 

1.197 

1.319 

1.711 

0.3 

1.052 

1.065 

1.082 

1-105 

1.141 

1.202 

1.328 

1.733 

0.4 

1.053 

1.066 

1.083 

1.107 

1.144 

1.206 

1.335 

1.750 

0.5 

1.054 

1.067 

1.084 

1.109 

1.146 

1.209 

1 .J 40 

1-764 

0.6 

1.054 

1.068 

1.085 

1.110 

1.148 

1-212 

1.344 

1.775 

0.7 

1.055 

1.068 

1.086 

1.111 

1,149 

. 1.214 

1.348 

1.785 

0.8 

1.055 

1.069 

1.087 

1.112 

1 . 151 

1.216 

1.351 

1.793 

0.9 

1.056 

1.069 

1.087 

1.113 

1.152 

1.217 

1.354 

1.800 

1.0 

1.056 

1.070 

1.088 

l: 113 

1.152 

1.219 

1.356 

1.806 
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TABLE 4-IV 

FACTOR FOR B = 0.030 AND F/2 = 0.970 


\ 

0.60 

0.65 

0. 70 

0.75 

0.80 

0.85 

0.90 

0.95 

0.1 

1.074 

1.092 

1.116 

1,150 

1.201 

1.289 

ll470 

2.065 

0.2 

1.077 

1.095 

1.120 

1.155 

1.209 

1-300 

1,491 

2.117 

0.3 

1.079 

1.098 

1.123 

1.159 

1.215 

1.309 

1.505 

2- 156 

0.4 

1.080 

1.099 

1.126 

1-162 

1.219 

1.315 

1.516 

2.185 

0.5 

1.081 

1.101 

1.127 

1.165 

1.222 

1.320 

1.525 

2.209 

0.6 

1.082 

1.102 

1.129 

1.167 

1.225 

1.324 

1.532 

2.229 

0.7 

1.083 

1.103 

1.130 

1.169 

1.227 

1-328 

1.538 

2.246 

0.8 

1.084 

1.104 

1.131 

1.170 

1.229 

1.331 

1.543 

2.261 

0.9 

1.084 

1.105 

1.132 

1.171 

1.231 

1.333 

1.548 

2.273 

1-0 

1.085 

1.105 

1 . 133 

1.172 

1.232 

1.335 

1.552 

2.285 
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TABLE 4-V 


FACTOR FOR B = 0.040 AND F/2 = 0.960 


\ ^ 
\to 

w \ 

0.60 

0.65 

0.70 

0.75 

0.80 

0.85 

0.90 

0.95 

0.1 

1.099 

i. 123 

1.155 

1.201 

1.271 

1.391 

1.641 

2.476 

0.2 

1.103 

1.128 

1.161 

1.209 

1.282 

1-407 

1.670 

2.555 

0.3 

1.105 

1.131 

1.166 

1.215 

1.290 

1.419 

1.691 

2.612 

0.4 

1.107 

1.133 

1.169 

1,219 

1.296 

1.428 

1.707 

2.657 

0.5 

1.109 

1.135 

1.171 

1-222 

1.300 

1.435 

1.720 

2.694 

0.6 

1.110 

1.137 

1.173 

1.225 

1.304 

1.441 

1.730 

2.725 

0.7 

1.111 

1. 138 

1.175 

1.227 

1.307 

1.445 

1.739 

2.751 

0.8 

1.112 

1.139 

1.176 

1.229 

1.310 

1.450 

1.747 

2.773, 

0.9 

1.113 

1.140 

1.178 

1.231 

1.312 

1.453 

1.753 

2,793 

1.0 

1.113 

1.141 

1.179 

1.232 

1.314 

1,456 

1.759 

2.810 
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TABLE 4-VI 

FACTOR FOR B = 0.050 AND F/2 = 0.950 


V 0/ 

\ to 

Nc 0 

0.60 

0.65 

0.70 

0.75 

0.80 

0.85 

0.90 

0.95 

0.1 

1.124 

1.154 

1.195 

1.253 

1.342 

1.495 

1.818 

2.914 

0.2 

1.129 

1.160 

1.203 

1.264 

1.357 

1-517 

1.857 

3.022 

0.3 

1.132 

1.165 

1.208 

1.271 

1.367 

1.532 

1.885 

3.101 

0.4 

1.135 

1.168 

1-213 

1.276 

1.374 

1.544 

1,906 

3. 164 

0.5 

1.137 

1.170 

1.216 

1.281 

1.380 

1,553 

1.923 

3.216 

0.6 

1.1-38 

1.172 

1.218 

1.284 

1.385 

1.561 

1.938 

3.259 

0.7 

1.140 

1.174 

1.221 

1.287 

1.389 

1.567 

1.950 

3.296 

0.8 

1.141 

1.176 

1.222 

1.290 

1.393 

1.573 

1-.960 

3.328 

0.9 

1.142 

1.177 

1.224 

1.292 

1.396 

1.578 

1.969 

3.357 

1.0 

1.143 

1.178 

1.226 

1-294 

1.399 

1.582 

1.977 

3,381 
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TABLE 4-VII 


FACTOR FOR B = 0.060 AND F/2 = 0.940 


\ W 

Xv 0 

P X 

0.60 

0.65 

0.70 

0.75 

0.80 

0.85 

0.90 

0.95 

0. 1 

1,150 

1.186 

1.23 6 

1.306 

1.415 

1.603 

2.002 

3.376 

0.2 

1.155 

1.194 

1.245 

1.319 

1.433 

1.630 

2.051 

3.517 

6.3 

1.159 

1.199 

1.252 

1.328 

1.445 

1.649 

2.087 

3.622 

0,4 

1.162 

1.203 

1.257 

1.335 

1.455 

1.664 

2.114 

3.705 

0.5 

1.165 

1.206 

1-261 

1.340 

1.462 

1.675 

2. 136 

3-774 

0.6- 

1.167 

1.208 

1.264 

1.344 

1.468 

1.685 

2. 155 

3.832 

0.7 

1.169 

1,210 

1.267 

1.348 

1.474 

1.693 

2. 171 

3.882 

0.8 

1.170 

1.212 

1.269 

1.351 

1.478 

1.701 

2. 184 

3.925 

0.9 

1.171 

1,214 

1.271 

1.354 

1.482 

1.707 

2. 196 

3.963 

1.0 

1.172 

1.215 

1-273 

1.356 

1.485 

1.712 

2.206 

3.997 
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It is apparant that for smaller values of 
(=0.80), the difference in the METREK and 
order of 15 percent. At higher values of 
the difference in the two results becomes 
order of many magnitudes. 


B(=0.02), as well as 
Gordon's results is of the 
.(=0-95), and B(=0.06)-, 
significant and is of the 


It must be noted that the value m =0.95 signifies almost a 

0 

totally scattering medium and in the absence of no absorption, one 
may expect to receive back all the radiance incident on the surface 
of the water. Such situations may occur in the middle of the ocean 


where the waters are almost turbidity free. For the LaRC experimental 
facility the values of larger than 0.80 are not expected. The 
water tank bottom at the LaRC's facility absorbs- some of the light 
photons, and remains -mostly imder 0.80, in which case the results 

f\j 

computed for > 0.80, and presented in this section may become 
meaningless. 


In conclusion, in this report: 

1) We have presented a matheioatical formalism applicable 
to the radiative transport process for the geometry 
of the light scattering relevant to the LaRC labora- 
tory experiment set up. Then an analytical solution 
to this equation has been derived to determine the 
backs cattered radiance as a function of the single 

'Vr 

scattering albedo o) , incident radiance I , and the 

0 o 

scattering phase function P(y), 

2) The results measured as well as calculated scattering 
phase functions for a wide range of turbid water 
conditions are summarized to determine the lower and 
upper bounding scattering phase functions. The 
bounding scattering phase functions have been used to 
obtain a better estimate of the ranges of the para- 
meters involved in the scattering phase function 
formulation. 


4-18 



ORIGINAi; PAGE IS 
OP POOR QUALITY 


3) The solution derived for the radiative transport 
equation is numerically evaluated for a range of 
scattering phase functions, and also for a range 
of single scattering albedo. 

Finally, 

4) METREK’ s results are compared with the already 
reported results in the literature on quasi-single 
scattering approximation. 

The differences between the two results have been discussed 
in Section 4.1. 
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APPENDIX A 


LOBATTO INTEGRATiON COMPUTER CODE 


A-1 



C HAIN PBOGEAH ‘QUADD* 

C THIS PEOGEAa PEEFOfiaS THE DESIEED INTEGEATION AND DETERMINES 

C the tESULTANT RADIANCE RECEIVED BY THE SEHOTE SENSOB- 

C 

C LOBATIO IMIEGBATION 
C DESCBIPIICN OF PABAMETEBS 
C A- LOWER LIMIT 

C B- DPPEE LIMIT 

C E- PRECISION 

C FCT- SUBROUTINE FCT(X,Y) 

c y=F (X) 

C USE EXTERNAL FCT IN MAIN PROGRAM 

C N- SET EQUAL TO 200 

C AUX- STORAGE FOE CALCULATED VALUES OF FCI,DJMENSXON TO 400 

C ICHEcK- <0, NOT CONVERGED WITH REQUIRED PRECISION 

C >0, CONVEIGED, NO. OF CYCLES REQUIRED 

C V- ANSWER 

C ARRAY- ARRAY OF PARAMETERS IN FCT, USE COMMON ARRAY IN FCT 

C IPARA- NO. OF PAEAHEIEBS IN ARRAY 

C 

IMPLICIT R£AL#3 (A-H,0-Z) 

DIMENSION AUX{5000) ,BBB (7) , TAD ( 10 , 8) ,AU(10) ,RAD(10,8) ,AAA(8) 
COMMON/XXX/ARRAY (20) 

PI=3. 14159 
BBB(1)=0.007 
BBB (2) =0.01 
BBB(3)=0.02 
BBB (4) =0.03 
BBB{5)=0.04 
BBB (6j=0.05 
BBB(7)=0.06 
AAA (1) =0.60 
AAA (2)=0.65 
AAA (3)^0. 70 
AAA (4) =0.75 
AAA (5) =0.80 
AAA(6)=0.85 
AAA (7) =0.90 
AAA(8)=0.95 
DO 35 J=1,7 
BB=BBB(J) 

F=2* (1.-BB) 

FF=F/2. 

WRITE (6,99) BB,FF 

99 FORMAI(*1' ,///, 15X, *B = ' , F6. 3 ,5X, ' F/2 =• ,R6.3,////) 

L=0 

DO 25 K=1,8 
L=L + 1 
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OMEGa=&AA (K) 

E=0. 00001 
DO 25 1=1,10 
IPAEA=4 

a=o. 000001 

B=PI/2. 

E=0. 000001 

ASSAY (1 ) =F10AT (IJ /10. 

ASSAY (2) =BB 
ARRAY 13>=0H£GA 
ARRAY (4) =F 
kU{l) =FLOAT (I)/10. 

EXTESMAl FCT 
H=2000 

CALL SPINT(A,B,E,FCO:,N,AOX,ICHSCK,V,HODE) 
AS=-ASSAY(1) =i'V/PI 
S1 =DEaP (AS) 

S11=S1* ((1.-0HEGA*E/2.) **1. 5) / ( (1 .-OMEG A) #*1.5) 
S2=1/l1. + ABRAY(l) ) 

S3=S2/ (1.-0MEGA*E/2.) 

S4=S3*0HEGA*BB/2. 

S5=S4*S11 
TAU (I,L)=S1 1 
SAD (I,L)=S5 

WRITE (7, 11 5) ARRAY (1) ,S5 

115 EOEHAI(2E8.5) 

25 CONTISOE 

DO 116 1=1,10 

116 WRITE (6, 117) AO (I) , (TA0(I,L) ,L=1,8) 

117 E0EHAT(2X,F5. 1,8 (2X,E8- 3) ) 

WRITE (6, 11 9) 

119 FOE HAT (////) 

DO 118 1=1,10 

118 HEITE(6,117) AU (I) , (SAD(I,L) ,L=1,8) 

35 CONTIROE 

STOP 

END 



no no on 


SOBEOUTIHE SPINI (A,B ,'E ,PCT,N , AOX, ICHECK, V, HOD'E) 
I.1P1ICIT AEAL*8 (A-H,0-Z) 

DIHENiilON A0X(1) 

1 IF {K.J.T. 12) N=12 
IF(A.EQ.B) GOTO 650 
IK= (N-4)/8 
IDV=2 
IX=IK+2 
IP=2*lK+3 
IFM=4*IK+3 
IFE=6=i=IK+3 
VC=1- iiO 
AA = A 
BB=B 

IF(A.LT.B) GOTO 10 

AA=B 

BD=A 

VC=-1-D0 

10 A0=0.U4761904761S0476D0 
A 1 = 0. ^7 682 6047361565900 
A2=0.43 174538120S8626D0 
A3=0- 487619047619047600 
T5=0. »30223 89627856 69D0 
T4=0. 468848793470714200 
T2=-i’4 
T1=-To 


FIRST INTERVAL (A TO B) 

C=0. 50* (BB-AA) 

D=AA+C 

X1=C*T1+D 

X2=C*T2+D 

X4=C*T4+0 

X5=C*‘x5+D 

CALL irCT(AA,AUX (IFE+1) ) 

CALL rCT(X1,F1) 

CALL FCT(X2,F2) 

CALL FCT (D,ADX (IFM+1) ) 

CALL FCT(X4,F4) 

CALL FCT(X5,F5) 

CALL rCT (BB, AUX (IFE+2) ) 

VO=C* lAO* (AUX (IFE+1) +AOX (IFE+2) ) +A1*(F1 +F5) +A2* (F2+F4) +A3*AUX (IFM+ 

11 )) 


SPLIT FIRST IMTBHVAL 
K=1 
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C=0. 5i)0#c 
V=0. DO 

AUX (IrE+3) =AUX (IFE+2) 

A'JX (lPE+2) =AUX (IFM+1) 

DO 100 1=1,2 
D=AA+FT,OAT (I+I-1) *C 
Xl = C*n+D 
X2=C*x2+D 
X4=C^i’4 + D 
X5=C^i5+D 
CALL £'CTtXl,F1) 

CALL fCT(X2,F2) 
lFHl=iFM+I 

CALL FCT {D,ADX (IFMI) ) 

CALL iCT(X4,F4) 

CALL iCT(X5,F5) 

IFEI=1FE+1 

AUX(Ii^l)=C*(AO+ (AUX (IFEI) +AUX (lFEl+1) ) +A1* (F1+F5) +A2*(F2+F4) +A3*A0 
1X(IF:ix)) 

100 Y=V+AUX (IPI) 

C 

C CHECK CONVERGENCE 

IF(V.iO,V0) GOTO 5000 
AUX (IOV+1) =V-V0 
VA=DABS (V) 

IF(VA.EQ.O.O) VA=DABS(V0> 

IF (DABS (AUX (IDV+1) } /VA.LT.E) GOTO 9000 
IF(IK.EQ.I) GOTO 500 
C END OP INITIAL CALCULATIONS 

C 

c BEGIN SELECTIVE SPLITTING PROCESS 

SDS=AOX (IDV+1) 

AUX(IX+1) =AA 
A0X(IX+2) =BB 
IT=1 
C 

C RESET VALUES 

1000 K=K+1 

IB=IT+1 

I1=IB+1 

IXIT=1X+IT 

C=0. 125D0* (AUX (IXIT+1) -AUX (IXIT) ) 

vo=v 

IDVIT=IDV+IT 

SDS=SDS-AUX(IDVIT) 

IH=IT+IT-1 

IPIR=1P+IR 
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V=V-AUX (IPIS) -AUX (IPIE+1) 
IXK=IX+K 

AaX(IXK+1)-=AUX (-IXK) 

KTW=2*K+1 

IF3KT«=IFE+KTW 

AUX (IFEKTH) =AOX (IFE^CTH-2) 

IF(K.EQ.IB) GOTO 105 

DO 104 I=I1,K 

L=K+I1-I 

1XL=IX+1 

AUX (IXL) =AUX (IXl- 1) 
IDVL=1DV+L 

AUX<IDVL) =AUX(IDV1-1) 

IZ=2*L-1 

IFEIZ=IFS+IZ 

AUX (liEIZ) =AUX (IFEIZ-2) 

AUX(IFEIZ+1) =AUX(IFEIZ-1) 

IFMIZ=IFM+IZ 

AUX =AUX (IFriIZ-2) 

AUX (Ii.1IZ+1) =AOX (IFMIZ-1) 
IPIZ=1?+IZ 

AUX <IPIZ) =AUX (IPIZ-2) 

104 AUX (IFIZ+1)=AOX(IPIZ-1) 

105 AUX (1X1T+1)=AUX (IXIT) +4.*C 

SPLIT LOOP 
IFEIS=1FS+IR 

AUX (IFEIB+2) =AUX (IFEIR+1) 

iFHia=iFa+iH 

AUX (IFSIS+ 1) =AUX {IFMIS) 

AUX (IiEIR+3) =AUX (IFXIR+I) 
AUX (1)=AUX (IPIR) 

AUX (2) --AUX (IPIH+1) 

L=0 

DO 300 J=1,2 
SA=0- DO 
DO 200 1=1,2 
L=L+1 

LV=IT+IT-2+L 

D=AaXtIXIT) +FLOAT (L+L-1) *C 

X1=C*T1+D 

X2=C*T2+D 

X4=C*I4+D 

X5=C*T5+D 

CALL FCTtX1,Fl) 

CALL FCT(X2,F2) 

IFMLV=IFB+LV 

CALL iCT(D,AUX (IFMIV) ) 
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CALL iCT{X4,F4) 

CALL FCT(X5,F5) 

IFELV=1FE+LV 

AT=C*iA0* (AUX (IFELV) +AUX (IFELV+1) ),-+Al* [FI +F5) +A2* (F2+F4) +A3*AUX (IF 
1MLV) ) 

SA=SA+AI 
IPLV=iF+LV 
200 A0X(IFLV)=AI 
V=V+.Sa 
JIT=-IT+ v7-1 
IDVJT‘i=IDV+Jir 
AUX (IDVJITJ =SA-AUX (J) 

300 SDS=Sj3S+AOX(1DVJ1T) 

CHECK CONVERGENCE- 
IF(V.EQ.VO) GOTO 9000 
VA=DAoS(V) 

IF{VA.EQ.0.0) VA=DABS(V0) 

IF(DAcS (SDS) /VA.LT.E) GOTO 9000 
1F(K, EQ.IK) GOTO 500 
C 

C FIND j.AEGE§T DELTA (WHERE TO SPLIT NEXT) 

DVHAX=DABS (AUX (IDV+1) ) 

DO 400 1=1, K 
IDVI=xDV+I 
ADV=DABS (AUX(IDVI) ) 

IF (ADV.LT. DVIIAX) GOTO 400 

DVHAX=ADV 

IT=I 

400 CONTINUE 
GOTO 1000 
9000 V=V*VC 

ICHECK=K 
RETURN 
500 V=V*VC 

ICHECK=-K 
RETURN 
650 V=0-D0 
RETURN 
END 
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SaBBOUTINE PCT-(X,-Y) 

IMPLICIT EEAL*8 (A-H,0-Z) 

COMMON/XXX/ABRAY (20) 

ZI^AEBAY (1) ARRAY (1) *DSIH (X) *DSIH (X) +DCOS (X) *DCOS (X) 
Z2=AERAY (2) *ARHAY (3) / (1 - -ARRAY (3) *ARRAY (4) /2. )' 
Z3=X/DTAN(X) 

Z4=DLOG (1.-Z2=®=Z3) 

Y=Z4/Z1 

RETURR 

END 



REFERENCES 


(1) A.H, Ghovanlou, J.N. Gupta, and R. G. Henderson, "Determination 
of Scattering Functions and Their Effects on Remote Sensing 
of Turbidity in Natural Waters,," The MITRE Corporation, NASA 
CR-145239, July 1977. 

(2.) A.H. Ghovanlou, J.N. Gupta, R.G. Henderson, and L. Poole, 
"Radiative Transfer Model For Remote Sensing of Suspended 
Sediments in Water," presented at the 4th Joint Conference 
on Sensing of Environmental Pollutants, November 6~11, 

New Orleans, LA. . . 

(3) A.H. Ghovanlou, et.al., "Laser Transmission Studies of East 
Coast Waters, Technical Report No. 2, sponsored jointly by 
Office of Nayal Research, Geography Branch, and The National 
Ocean Survey, Engineering Department Laboratory, and The 
U.S. Geological .Survey, March 1973. 

(4) ,A.H. Ghovanlou, "Analytical Model For Remote Sensing of 
Water Turbidity," The MITRE Corporation, NASA CR-145050, 
September 1976. 

(5) H.R. Gordon, "Simple Calculation -of the- Diffuse Reflectance 
of The Ocean," Applied Optics, Vol. 12, No. 12, December 1973. 

(6) V. Granatstein, M. Rhinewine, A. Levine, D. -Feinstein, M. 
Mazurowski, and K. Piech, "Multiple Scattering of Laser 
Light Fram A Turbid Medium," Applied Optics, Vol. 11, 

No. 5, May 1972. 

(7) K.M. Watson, J. Math, Phys., 10, 688, 1969. 

(8) y.N. Barabanekov and V. Finkelbery, Soviet Physics JETP 26, 

587, 1968. 

(9) G. MiTtonwing, An Introduction To Transport Theory, John 
Wiley and Sons, Dec., New York, 1962. 

10) E. Titchmarsh, The Theory of Functions,, The Clarendon Press, 
Oxford, 1932. 

11) R.V. Churchill, Complex Variables and Applications, McGraw 
Hill Book Company, New York, 1960. 

12) S.Q. Duntly, "Light in the Sea," J. Opt. Soc. America, Vol. 53, 
No. 2, February 1963. 


B-1 



REFERENCES (Concluded) 


(13) G- Kullenberg, Rep 'Inst, Fys , Oceanog. , 5, 16. Kobehhauns 
Universitet. , 1969. 

(14) N.G, Jerlov, Medd. Oceanog, Inst. Goteborg, 30, 1961. 

(15) J-E. Tyler, Liiran, and Oceanog., 6, 1961. 

(16) G. Kullenberg and 0. Berg, Rep. Inst. Fys. Oceanog., 19, 40, 
Kobenhavns Universitet, 1972. 

(17) G. Kullenberg, Deep Sea Res., 15, 1968. 

(18) T.J. Petzold, "Volume Scattering Functions for Selected Ocean 
Waters," AD-753474, Scripps Institution of Oceanography, 

October 1972. ' • 

(19) M. Abramowitz, and I. Stegum, Editors, "Handbook of Mathematical 
Functions," Dover Publications, Inc., New York, 1970. 


B-2 



DISTRIBUTION LIST 


A-10 

G. MacDonald 

D-10 

■S . Blum 
W. Gouse 
C. Zraket 

D-11 

W. Sievers 
D-12 

H. Benin gt on 

C. Grandy 

A. Tachmindji 

W-50 

R. Greeley 
R. Ouellette 
J. Golden 
R. Pikul 

R. Foreman 

D. Sluyter 

W-51 

E. Sharp 

L. Gsellman 

W-52 

M. Scholl 

G, Bennington 

S . Goldstein 


W-53 

T. Wright 
E. Ward 

Ali Ghovanlou (10) 

W. Flury 
T. Kuch 

M. Barbier 
0 . Farah 

E. Friedman 
J. Gupta (10) 

E. Keitz 

R. Henderson 

N. Lord 

W-54 

E. Krajeski 
W-55 

A. Chain s** 

N. Coates 
J. Stone 

W-56 

L. Thomas 

S. Lubore 
G. Erskine 

W-57 
S . Lewis 

W-50 Library (5) 

Metrek Library 

Technical Report Center (2) 


C-1 


Document Control (8) 
Westpart Storage* (15) 



EXTERNAL DISTRIBUTION LIST 

No. of 
Copies 

NASA Langley Research Center 
Hampton, VA 23665 

Attn: Report & Manuscript Control Office 

Mail Stop 180A 

Lamont R. Poole, Mail Stop 272 

2 

59 

NASA Ames Research Center 
Moffett Field, CA 94035 
Attn: Library, Mail Stop 202-3 

1 

NASA Dryden Flight Research Center 
P. 0. Box 273 
Edwards, CA 93523 
Attn: Library 

1 

NASA Lyndon B. Johnson Space Center 
2101 Webster Seabrook Road 
Houston, TX 77058 
Attn: JM6 /Library 

1 

NASA Marshall Space Flight Center 
Marshall Space Flight Center, AL 35812 
Attn: Library, AS61L 

1. 

•Jet Propulsion Laboratory 

4800 Oak Grove Drive 

Pasadena, CA 91103 

Attn: Library, Mail Stop 111-113 

1 

NASA Lewis Research Center 
21000 Brookpark Road 
Cleveland, OH 44135 
Attn: Library, Mail Stop 60-3 

1 

NASA John F. Kennedy Space Center 
Kennedy Space Center, FL 32899 
Attn: Library, NWSI-D 

1 

National Aeronautics & Space Administration 
Washington, D.C. 20546 
Attn: ER 

1 

C-2 




EXTERNAL DISTRIBUTION EIST (Concluded) 


No. of 
Copies 


NASA Scientific & Technical Information Facility , 
6571 Elkridge Landing Road 
Linthicum, Heights, MD 21090 


30 



